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We consider a two-component Fermi gas interacting via a Feshbach molecular state. It is shown 
that an important energy scale is Eg = g*m'^/(647r^) where g is the Feshbach coupling constant 
and m the mass of the particles. Only when Eg ^ ep where ep is the Fermi energy can the gas be 
expected to enter a universal state in the unitarity limit on the atomic side of the resonance where 
there are no molecules present. The universal state is distinct from the molecular gas state on the 
other side of the resonance. We furthermore calculate the energy of the gas for this universal state 
and our results are related to current experiments on ®Li and ^''K. 



] One of the most interesting recent developments within 
the field of ultracold atom gases is the use of Feshbach 
resonances to manipulate the effective atom-atom inter- 
action. By controlling an external magnetic field, the 
Feshbach molecular state can be tuned to an energy close 
to the scattering atoms. In this case, the vacuum atom- 
atom scattering is unitarity limited with a diverging scat- 
tering length a and the system is strongly interacting. 
Thus, it is possible with the atomic gases to address the 
important question of the nature of a many-body system 
with an interaction characterized by \a\ —> oo. Indeed, 
remarkable experimental progress concerning this ques- 
tion has been made in recent months with several groups 
reporting results for two-component atomic Fermi (^Li 

, and ''^"K) gases interacting via a closed channel Feshbach 
molecular state close to resonance Q, |^ IE Q - 

In this paper, we examine the problem of a two com- 
ponent Fermi gas interacting via a Feshbach molecular 

I state close to zero energy. We show that only for suffi- 
ciently broad resonances with Eg 3> cf can the gas be 

[ expected to exhibit universal behavior close to resonance 
on the a < when there are no molecules. The possi- 
bility of such universal behavior when ]a| — > oo has been 
examined by a number of authors [E B 13 ■ The ba- 
sic idea is that \a\ is not a relevant length scale in this 
limit and that the thermodynamic quantities of the gas 
therefore should only depend on the density n and the 
temperature T. The gas is in this sense universal since 
its properties are independent of the specific details of 
the interaction. The energy per particle for a Fermi gas 
with T <^T-p where T-p is the Fermi temperature is pre- 
dicted to scale as E/N — 3eF(l + f3)/b in the universal 
limit and various values for the constant (3 has been pre- 
dicted 0, S Hi- For a Bose gas with \a\ oo, it was 

\ argued that E/N oc 0- Similar results where de- 

rived for bosonic three particle states ■ We will in the 

' present paper calculate (3 for a Fermi gas in the universal 
limit for Eg ^ ep taking into account two-body scatter- 
ing in a medium. Our results are finally related to the 
present experiments on the atomic gases. 

The explicit model considered here consists of atoms 
with mass m in two hyperfine states denoted ] and [ be- 
ing coupled to a closed channel with a Feshbach molecule. 
The Hamiltonian describing this interacting mixture of 



atoms and molecules (described as point bosons) is 

K,q ^ ''^ 

+ 2^ y «k+qT"k-qi«k-q'iak+q'T 

k,q,q' 

with iJo = Ek,.efcaLaka + Ek^k'^k^k. Here 
^imic _ 2vha,vc + K'^/4:m, where 2z^barc is the energy of 
a bare molecule with momentum zero measured with re- 
spect to the energy of a pair of atoms at rest in the open 
channel, e/c = k"^ /2m is the kinetic energy of an atom, 
gharcio) is the bare molecule-atom coupling matrix ele- 
ment, V is the non-resonant interaction between atoms, 
and V is the volume of the system. 

The main object to calculate here is the thermody- 
namic potential for the gas = —ksTlnZ with Z = 
Tr{exp[-/3(H-AiNF-2^NB)]}. Here, Np = Ek,>L«k. 

and Nb = X^k ^k^k the number of bare atoms and 
molecules respectively and /i is the chemical potential. 
We assume an equal density of the two atomic hyper- 
fine states and a total density n ~ fcp/37r^. Also, the 
molecules and atoms are assumed to be in equilibrium. 
Using the ladder approximation (or the Gaussian approx- 
imation for the partition function), the thermodynamic 
potential can after some algebra be expressed as 

VI r (PK r duj 1 

V ^ V (27r)3 J ~ e'^'^ -1^ 

ImTrln[l - (V + gbarcDo(K, u;)gbarc)G[,') (K, w)]. (1) 

Here fio is the thermodynamic potential of an ideal 
gas of non-interacting atoms and molecules described by 
Hq. We omit for brevity here and in the rest of the 
paper a small positive imaginary part to the frequen- 
cies in the (retarded) Greens functions which should re- 
ally be evaluated aX oj + i6. The free molecule Greens 
function is Dq{K,uj)^^ = lo — £'^^™ -|- 2/i and the two 

particle Greens function Gq (K, q, cj) — [/(CK/2+q) + 
/(fK/2-q) - l]/KK/2+q + fK/2-q " dcscribcs thc free 
propagation of a pair of atoms with COM momentum K 
and relative momentum q in the presence of the medium. 
Here = Ck — fJ. and f{x) = [exp{x/kBT) + is the 
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Fermi function. The trace in Eq. refers to integration 
over the momentum variables appearing in the matri- 

ces F(q,q') and gbaro(g)G^ ffbare(q'')- Various forms of 
Eq. ^ have been discussed in the Hterature in connec- 
tion with the BCS-BEC crossover problem For the 
present purposes, it is convenient to rewrite Eq. |^ as 
17 = Jlo + AfJ with 



An 



(2^ 



duj 



{ImTrln[G[,'^(K,c.)/G^^g^(K,c.)] 



(2), 



+ImTrln[Do(K,cj)/D(K),w)]}(e'^/^«T _ (2) 

where Gf^g = Gq —V describes the propagation of 
a pair of atoms interacting only via the background in- 
teraction in the ladder approximation. D^^ = D^^ — 11 
is the molecule Greens function with the self energy 

Ii{K,Lu) = Tr[gbaroG[,|^(K,L^)gbare] [H ■ Equation © is 
convenient since it separates the corrections to coming 
from the Feshbach molecule from the effects of the back- 
ground (non-resonant) interaction. Since we focus in the 
effects of the resonant part of the interaction, we will in 
the rest of this paper only consider the molecule term in 
Eq. 0. The \tl[Dq{K,uj)] term in Eq. Q is trivial as it 
simply cancels the bare molecule contribution to rig- 

In order to be able to express our result in terms of 
observables only, we use the low energy effective theory 
developed in ref. We therefore write zD~^{K,lj) — 
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K 



ig rrv 



3/2 



(47r)-iA 



— 2fi and uj 



- Il{K,u) with Ek = 
2fi - K^/im. Here 2^ is 
the energy of the molecule taking into account the high 
energy dressing by the atoms but excluding the threshold 
effects which are explicitly treated by the term. The 
wave function renormalization of the molecule given by 
z is due to the coupling to high energy atoms. The con- 
stant g yields the coupling between this dressed molecule 
and the atoms and 11 gives the medium corrections to 
the molecule self energy 11. Explicit expressions for these 
quantities are given in ref. |12i | and will not be repeated 
here. The point is that g and 21/ can be extracted from 
experiments without any arbitrary assumptions of the 
behavior of the bare quantities. 

It is easy to show that Eq. Q reproduces the two- 
body term in the virial expansion of in the Boltzmann 
regime |l3j . This expansion has been used to prove that 
a gas of atoms interacting via a closed channel resonant 
(zero energy) Feshbach molecule exhibits universal be- 
havior in the non-degenerate limit [q]. We will presently 
examine under which conditions such universal behavior 
persists in the low temperature degenerate regime. 

We split Af2 into a contributions coming from un- 
damped and damped two particle states. This is done 
by writing Eq. Q (neglecting the background parts) as 
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arctan 



d^Kd^q /(gK/2+q)/(eK/2-q) 
(27r)6 Imn(K, eK/2+q + U/2-^)/g^ 

Imn(K, ^K/2+q + $K/2-q) 
g2/m -2l/- Rcn(K, CK/2+q + ^K/2-q) 



where A^m = keT J2k'^og[l - exp(c^/f /fc^T)] [1 - 
dw^{k,LL!K)] is the contribution from the undamped 
molecular states with energy luk- The energies luk = 
ujk + K'^/^m — 2yU are determined by ujk ^ 2v + 
g'^m^/'^('iTT)-^^f^u)]^+ii{K,u)K) withtj^ < K'^/Am~2^ 
since Imn(K,[x)) 7^ for u; > /Am — 2/i. It is easy to 
show that IV{K,uj) > for w < K^/Am - 2/Lt, i.e. the 
medium only increases the molecule energies. There will 
therefore be no undamped molecules below the Fermi sea 
for v > 0. To arrive at the second term in Eq. from 
Eq. 121 , we have divided and multiplied by Imll and per- 
formed the Lo integration. Using Eq. (^J, we can analyze 
under which conditions the interacting gas will exhibit 
universal thermodynamics. 

As explained above, the universal thermodynamics hy- 
pothesis states that the thermodynamic quantities should 
depend only on T and n in the unitarity limit when 
\a\ — > 00. For the specific system considered here, the 
scattering length is ATraics/m = —g^/2i'. The thermody- 
namic potential depends therefore in general on fijV^T, 
and g and we write n(T,fj,, i^,g)- If such a gas is to ex- 
hibit universal thermodynamic behavior, we must have 
n{T, g) ~ r2(T, /i) close to resonance, i.e. the ther- 
modynamics become independent of the specific details 
of the resonance. 

Let us first briefly consider the molecular side of the 
resonance with v < 0. Since fJ. < i^^ the system in 
this case consist of a pure molecular gas with essen- 
tially no atoms present for T <^ Tp and we need only 
to consider Ai7j\/ in Eq. We therefore also have 

n = and the molecule energies are given by w = 
— {yjEg — 2v— ^JWgf' . We have introduced the energy 
Eg = g^w? I {(SA-K^) which determines the range of ener- 
gies for which threshold effects are important. So univer- 
sality for < is straightforward to understand within 
the present model: For —Eg ^ 2z^ < 0, we simply obtain 
an ideal gas of molecules with energy d) = —Ti'/{ma^^g). 
However, this result ignores the molecule-molecule inter- 
actions, which are not included in the present formalism. 
Such effects could be important for v < Q where there is 
a macroscopic number of molecules [T^ . 

The main purpose of the present paper is to examine 
the arcs < (i^ > 0) side of the resonance where there are 
no undamped molecules and AnM in Eq. Q can be ig- 
nored. The only place where g and v enter in the expres- 
sion for n is then in the argument of the arctan(. . .) where 
V appear explicitly and g^ enters since H oc g^. The 
molecule self energy H scales as ~ 0{g'^mkY / An) 12]. 
From this it follows that in order to make the argument 
of arctan(. . .) independent of g^ and i/, we must have 
g^rnkF/Air 3> ep- Only in this case will the H terms 
in the arctan(. . .) dominate and thus cancel the g^ and 
v dependence even when 3> ep- The reason for the 
requirement v ^ ep is that fj- < v from particle conser- 
vation. When 1/ < ep we have ^ ~ due to the presence 
of damped molecules; i.e. n and the thermodynamics of 
the gas will depend critically on v when < v ep- It 
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should be noted that for v/e^ 0, we of course approach 
the pure molecular gas state for i/ < described above; 
the state of the gas changes continuously across the res- 
onance. What we have shown here however is that the 
gas approaches a different universal state independent 
of (g, ly) for I' 3> ep only if the resonance is sufficiently 
broad with ^ AnkFrn^^- This state is characterized 
by strongly correlated atom pairs with no real molecules 
present. The condition for universality on the Oics < 
side can also be written 

Eg » ep. (4) 

In terms of the scattering length 4:TTa,.cs/m = —g^/2i^, Eq. 
(@J) can be written as fepjarcsl ^ 1 even when v ^ ep. 
Thus, the condition for universality on the atomic side 
of the resonance simply states that the coupling must 
be so strong that one can reach the unitarity regime 
fcplaresl ^ 1 before there is a significant population of 
molecules. The in the denominator of the arctan(. . .) 
may be regarded as arising from an effective range term. 
We can then identify the effective range for the effec- 
tive interaction mediated by the Feshbach molecule as 
Teff = —Sirm'^^g"^. Using this, we can finally write Eq. 
(0J) in the illuminating way /cp|rcg| ^ 1, i.e. universality 
can only be reached if the effective range of the interac- 
tion is small. 

We now consider the arcs < and T = limit of the 
universal behavior for in more detail. For T = 0, we 
can write Afi/V ~ ^n\i{p)I where liv j [i^g^rr?!'^ j ^fji) 
is the T = integral in Eq. (O expressed in dimcn- 
sionless variables (fc — k/ ^2m^ etc.) and nid(M) = 
(2m/x)^/^(37r^)~^. The integral / can be evaluated nu- 
merically. In the universal limit e^, we get /un; ~ 
-0.3 for T = 0. This immediately gives Vl{T = 0,/i) = 
-2^nid(//)(l - 5/uni/2)/5 and ^ = {2my^{iT:^nf'^{l - 
5/uni/2)"2/^. We finally obtain E/N = 3eF(l+/3)/5 with 

/3 = (1 - 54„i/2)-2/3 _ 1 ^ _o.3. (5) 

Equation Q gives the T = energy correction in the uni- 
versal limit when all two-particle correlations are evalu- 
ated using the ladder approximation. It should be noted 
that the medium effects giving energy shifts and Pauli 
blocking through 11 are significant; if the vacuum values 
Imllvac — g^mq/47r and 11 = are used in Eq. (|3Jl, we 
obtain /uni vac — 5/8 which gives (3 ~ —0.519. 

To illustrate the conclusions described above, we 
present results based on a numerical evaluation of Eq. 
0. In Fig. n we plot the energy per particle E/N as a 
function of — fcpaics for two different coupling strengths: 
^siargc > <^F and < ep. For giarge, we also 

plot the energy for two different temperatures T = 
and T — Tp/A. The number of particles is fixed while 
Ores = —g'^m/{%'Kv) is varied through varying v. We see 
that for (/large, the energy per particle approaches a con- 
stant for — A;pares ^ 1. In agreement with the Eq. (|3J), we 
find E/N ~ 0.42ep for -fepCros > 1 and T = 0. For the 
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FIG. 1: The energy per particle E/N in units of ep as a 
function of — fcrires for a narrow and a broad resonance. 

narrow resonance Ismail on the other hand, the energy per 
particle steadily decreases with increasing — fcparcs and 
the gas does not reach the universal limit for arcs < 0. 
This is simply because when — fcparos ?i 1, we already 
have V < ey for such a small value of g, and the de- 
pletion of the Fermi sea into the molecular level is sig- 
nificant. Of course, the depletion of the Fermi sea for 

< ep also happens for giargc- However, since this corre- 
sponds to — /cparcs ^ 5'fargo"^^/4'''^F' there is a large re- 
gion of parameter space 1 <C — fepOros ^ 5iaigc'™^/4'''fcF 
where the thermodynamics of the gas is universal with no 
molecules present. This is illustrated by the inset, which 
plots E/N for giarge over a larger range of — fepOres- In a 
second inset, we plot E/N for — fepCres 1- The solid 
line is the numerical results whereas the dashed line is the 
Galitskii result E/N = eF[3/5 + 2fcFarcs/(37r) + 4(11 - 
21og2)/(357r2)(A;Fa)2 + ...]. We see that our theory re- 
produces the Galitskii result for — fcparcs <C 1 as it should. 

The depletion effect for v e-p \s illustrated more 
clearly in Fig. [3 where we plot the chemical poten- 
tial /i and E/N for the same coupling strengths but 
now as function of v for T = Tp/A. We also plot the 
density of (damped) molecules rimoi = TrZ3 and atoms 
'T-atom = n — rimoi for giargc- Wc sce that there is only a 
significant number of molecules when v ep as expected. 
They will eventually become the undamped molecule gas 
for V <Q. Of course, for both coupling strengths the gas 
converges toward a pure molecular state for v <0. 

In total, the numerical results confirm our conclusion, 
that the thermodynamics of the gas only exhibit univer- 
sal behavior when — arcs^p 3> 1 for a sufficiently broad 
resonance with Eg^_^^ ^ ^ ep. 

Let us now compare our result for the T = universal 
limit with a number of recent calculations. In ref. 0, 
the value j3 ~ —0.67 was estimated from an expansion of 
the on-shell scattering amplitude in powers of fcpa and 
an approximate momentum average. In ref. 01, the val- 
ues (3 ~ —0.43 and (3 ~ —0.67 were suggested from two 
different Fade approximants of the ladder series for the 
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FIG. 2: The chemical potential /i, the energy per particle 
E/N in units of ep and the density of atoms and molecules in 
units of Tier = n/2 = (2meF)^^'^/(67r^) as a function of v for a 
narrow and a broad resonance. 

energy. The T = result j3 ~ —0.3 presented in this 
paper however, is the exact within the ladder approxi- 
mation excluding superfluid correlations. A fixed node 
Greens function Monte Carlo calculation with ^ 40 par- 
ticles yields (i ~ —0.46 and [3 ~ —0.56 when superfluid 
correlations are excluded and included respectively @. 
It could be interesting to examine in more detail the rea- 
sons for the difference between the present diagrammatic 
result and the Monte Carlo results in ref. 8]. One possi- 
ble reason could be > 2 body correlations which cannot 
be expected to be small in the unitarity limit. 

We conclude this paper by relating our results to some 
of the experiments on atomic gases. First, we exam- 
ine whether a couple of the experimentally relevant res- 
onances can be expected to exhibit universal thermody- 
namics. For one such Feshbach resonance at the magnetic 
field Bq ~ 200 G for ""^K in the two open channel hyper- 



fine states |9/2,-9/2) and |9/2,-7/2) [Ti|, an analysis 

— 1/3 

yields feplaresl — 15nn ^y/v where ni2 is the density in 
units of 10^^ cm^"^ ^M,- Another case of interest is the 
^Li Feshbach resonance between the two lowest hyperfine 
states at Bq ~ 860 G . An estimate of this resonance 

yields feplarcsl — 1-1 x lQ^rh[2^^ / Vhaxe E3- This shows 
that these two particular resonances are sufficiently broad 
(large g"^) such that one would expect universal behav- 
ior on the Ores < side of the resonance. There has 
recently been a number of measurements of the energy 
of a ^Li gas close to the Feshbach resonance considered 
above PjHj Q • In agreement with the conclusions of this 
paper, they suggest that the properties of the gas indeed 
become universal in the limit \a\ — > oo. The Duke group 
obtain (3 = -0.26 ± 0.07 on the a^s < side of the ^Li 
resonance llj . This value agrees well with Eq. |SJ| . How- 
ever, the Innsbruck group measures p = — 0.681q for 
the same resonance Q . Finally, the group at ENS obtains 
{Etot - i^kin)/i^kin ^ -0.3 for arcs -oo (3- But since 
both the total energy i^tot and the kinetic energy i?kin 
of the gas are measured for the interacting system this 
ratio is not identical to the parameter [3 as defined in the 
present paper. In general, further investigation is needed 
in order to remove some of the apparent inconsistencies 
of the results for \a\ — s- oo reported so far. 

In conclusion, we have shown that a two-component 
Fermi gas interacting via a Feshbach state can be 
expected to exhibit universal behavior on the Ores < 
side only when the resonance is sufficiently broad with 
Eg ^ ep. We furthermore calculated the universal 
parameter (3 taking into account two-body processes in 
the ladder approximation. Our results were compared 
to the experimental results when appropriate. 

We acknowledge very useful discussions with C. J. 
Pethick, H. Heiselberg, and S. J. J. M. F. Kokkclmans. 
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